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Abstract. This paper is concerned with the study of a geometric flow whose law involves 

... a singular integral operator. This operator is used to define a non-local mean curvature 

^j^ I of a set. Moreover the associated flow appears in two important applications: dislocation 

^^ ■ dynamics and phasefield theory for fractional reaction-diffusion equations. It is defined by 

fvj I using the level set method. The main results of this paper are: on one hand, the proper 

. ■ level set formulation of the geometric flow; on the other hand, stability and comparison 

C^ ' results for the geometric equation associated with the flow. 
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1 Introduction 



^ ! In this paper, we define a geometric flow whose law is non-local. We recall that a 

^ I geometric flow of a set fi is a family {flt}t>o such that the velocity of a point x G dflt 

^ I along its outer normal n{x) is a given function of x and n{x) for instance. In our case, this 

^ ■ velocity does not only depend on x and n{x) but also on a fractional mean curvature at x. 

Our motivation comes from two different problems: dislocation dynamics and phasefield 

theory for fractional reaction-diffusion equations. 



1.1 Motivation and existing results 



S^ . Mathematical study of non-local moving fronts recently attracted a lot of attention 

(see in particular [12] and references therein). An important application is the study of 
dislocation dynamics [3]. 

Dislocation dynamics 

Dislocations are linear defects in crystals and the study of their motion gives rise to 
the study of a non-local geometric flow. In recent years, several papers were dedicated to 
this problem. We next briefly recall the results contained in these papers. 

A dislocation creates an elastic field in the whole space M^ and this field creates a 
force (called the Peach-Koehler force) that acts not only on the dislocation that created 
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it (self-force) but also on all dislocations in the material. We restrict ourselves here to 
the case of a single curve. We also assume that this curve moves in a plane (called the 
slip plane). 

The level set approach [2U [131 tlH] is a general method for constructing moving inter- 
faces. It consists in representing Qt as zero level sets of functions u{t, ■). The geometric 
law satisfied by the interface dQt is thus translated into an evolution equation satisfied 
by u. This approach is used in [Sj to describe the dynamics of a dislocation line. If dQt is 
the zero level set of a function u{t,-), the following non-local eikonal equation is obtained 

dfU = {ci{x) + k[x, u])\Du\ 

where Ci is an external force and k[x, m] is the Peach-Koehler force applied to the curve 
(A^ = 2 in this application). 

We briefly mentioned above that the Peach-Koehler force is created by the curve. Let 
us be a bit more specific. This force is computed through the resolution of an elliptic 
equation on a half space (corresponding to the law of linear elasticity). This equation 
is supplemented with Dirichlet boundary conditions. On one hand, the boundary datum 
equals the indicator function of the interior of the curve. On the other hand, loosely 
speaking, the force on the curve equals the normal derivative of the solution of the elliptic 
equation. Hence, the integral operator which defines the Peach-Koehler force is a Dirichet- 
to-Neumann operator associated with an elliptic equation. In particular, the operator is 
singular. 

In order to define solutions for small times, the authors of [3] consider a physically 
relevant regularized problem and k[x, m] reduces to 

/ co{z)dz 

J{z:u{z)>0} 

with Co G 1^1'^ (M^). 

The major technical difficulty of this paper is that Cq does not have a constant sign and 

consequently, solutions corresponding to ordered initial data are not ordered; in other 

words, comparison principle does not hold true. In particular, this is one of the reasons 

why solutions are constructed for small times. If Ci is assumed to be large enough, Alvarez, 

Cardaliaguet and Monneau [2] managed to prove the existence and uniqueness for large 

times. 

The difficulty related to comparison principle is circumvented in [19] by assuming that 
the negative part of cq is concentrated at the origin. The Peach-Koehler force k[x, u] (in 
the case of a single dislocation line) is defined in [19] as 



sign(M(x + z) — u{x))cQ{z)dz = / co{z)dz — / co{z)dz (1) 

J {z:u{x+z)>u{x)} J {z:u(x+z)<u(x)} 

where sign(r) equals 1 if r > and — 1 if r < 0. After an approximation procedure, the 
problem can be reduced to the study of 



dtU 



ci(x) + / (f/(x + z)- U{x))co{z)dz 



\DU\ 



where Cq is smooth, non-negative and of finite mass. We used the letter U instead of u 
in order to emphasize the fact that a change of unknown function is needed in order to 
reduce the study of the original equation to the study of this new one. 

A second important remark is that solving such non-local eikonal equations does not permit 
to construct properly a geometric flow. More precisely, if the initial front dQo is described 
with two different initial functions Uq and Vq, it is not sure that the zero level sets of the 
corresponding solutions u and v coincide. In other words, the invariance principle does 
not hold true. 

Still assuming that the negative part of cq is concentrated at the origin, a good geo- 
metric definition of the fiow is obtained in [17] by considering a formulation "a la Slepcev" 
of the geometric fiow. The equation now becomes 



dtu 



ci(x) + 



co{z)dz 



{z:u{t,x+z)>u{t,x)} 



\Du\ 



(2) 



We point out that, with such a formulation, we cannot deal with singular potentials cq. 

Notice that in [17], several fronts move, and they are interacting. The motion of a 
single front is a special case. Eventually, existence results of very weak solutions in a very 
general setting are obtained in [S] and uniqueness is studied in [Hj. 



In [T3], it is proved that if Co{z) is smooth and regular near the origin and behaves 
exactly like |2;|~^~^ at infinity, then a proper rescaling of (E]) converges towards the mean 
curvature motion (see Proposition [1] and Theorem H] below). 

We finally mention that Caffarelli and Souganidis [H] consider a Bence-Merriman- 
Osher scheme with kernels associated with the fractional heat equation (that is to say the 
heat equation where the usual Laplacian is replaced with the fractional one). They prove 
that this scheme approximates the geometric flow at stake in this paper. 



Phasefield theory for fractional reaction-diflFusion equations 

Our second main motivation comes from phasefleld theory for fractional reaction- 
diffusion equations [21]. If one considers for instance stochastic Ising models with Kac 
potentials with very slow decay at inflnity (like a power law with proper exponent), then 
the study of the resulting mean fleld equation (after proper rescaling) is closely related to 
phasefleld theory for fractional reaction-diffusion equations such as 



dtu' + (-A)°/2^^ + 



rl+a 



f{u^ 







where (— A)"/^ denotes the fractional Laplacian with a E (0, 1) (in the case presented 
here) and / is a bistable non-linearity. In particular, it is essential in the analysis to deal 
with singular potentials. Indeed, we have to be able to treat the case where 



co{z) 



1 

\N+a 



with a G (0, 1). It is also convenient to use the notion of generahzed flows introduced by 
Barles and Souganidis [H] in order to develop a phasefield theory for such reaction-diffusion 
equations. See [21] for further details and [18] for analogous problems. 



1.2 A new formulation 

The main contributions of this paper are the following: 

• to give a proper level set formulation of dislocation dynamics for singular interac- 
tion potentials; in particular, sufficient conditions on the singularity to get stability 
results and comparison principles are exhibited; 

• to shead light on the fact that the integral operator measures in a non-local way the 
curvature of the interface; 

• to study the geometric flow in detail: consistency of the definition, equivalent defi- 
nition in terms of generalized fiows, motion of bounded sets etc. 

Because v{dz) = co{z)dz is singular, we cannot define k[x,m] as in ([2]). Indeed, we 
must compensate the singularity as it is commonly done in order to get a proper integral 
representation of the fractional Laplacian. We recall that the fractional Laplacian can be 
defined as follows 

(-A)"/2^(a;) = -cjv(a) [ {u{x + z) - u{x)) "^"^ 



^\N+a 



where CAr(a) is a given positive constant depending on N,a. Notice that if a < 1 and 
u is Lipschitz continuous at x and u is globally bounded, the integral is well defined. If 
a > 1, the integral is not convergent in the neighbourhood oi z = 0. In this case, the 
fractional Laplacian is defined either by considering the principal value of the previous 
singular integral or by writing 

(-A)°/2^(x) = -CN{a) / {u{x + z)- u{x) - Du{x) ■ z1b{z)) 



^\N+a 



where 1b{z) denotes the indicator function of the unit ball B. Notice that we used the 
fact that the singular measure 

^(dz) = —^ (3) 

(with a G (0, 2)) is even in order to get (at least formally) 

dz 
{Du{x)-z1b{z))-^ = 0. 

As far as the fractional mean curvature is concerned, we must compensate the singular- 
ity of the measure i/ in a geometrical way. We explain how to do it when ui^dz) = co{z)dz 



with co(^;) = |2;| ^ °. Hence, we start from ([T]). We use the fact that cq is even in order 
to get (formally) 

u{z e M^ : Du{x) ■ z>0) = u{z eR^ : Du{x) ■ z < 0) . 

Straightforward computations yield 



/ 



sign(ti(a; + z) — u{x))cQ{z)dz = 

^{z : u{x + z) > u{x), Du{x) ■ 2; < 0} — u{z : u{x + z) < u{x), Du{x) ■ z > 0} 



We thus define an integral operator /i;[x, n] for a general singular non- negative measure u 
as follows 

k[x, u] = u{z : u{x + z) > u{x), Du{x)-z < 0} — h'{z : u{x + z) < u{x), Du{x)-z > 0} . (4) 

We explain below in detail (see Lemma [2]) the rigourous links between the different for- 
mulations we considered up to now. 

Notice that this definition makes sense even if u is not even. We recall that the 
fractional Laplacian is a Levy operator. Since Levy operators |1] are defined for singular 
(Levy) measure that are not necessarily even, this seems to be relevant to define the 
fractional mean curvature for singular measures that are not necessarily even. 

We can say that this singular integral operator measures in a non-local way the curva- 
ture of the "curve" {u = u{x)}. Indeed, loosely speaking, we can say that in Formula (jlj) 
the first part (resp. the second one) measures how concave (resp. convex) is the set 
Q = {z : u{x + z) > u{x)} "near x". Moreover, we prove (see Proposition [2] below) that, 
when i/ is given by ([3]), the function {l — a)K[x, u] converges as a G (0, 1) goes to 1 towards 
the classical mean curvature of {u = u{x)} at x. This is the reason why we refer to k,[x, u] 
as the fractional mean curvature of the curve {u = u{x)} at point x. 

The variational case 

When the singular measure uldz) has the form 

u{dz) = -(v-G{zyjdz 

for a vectorfield G, the previous singular integral operator can be written as follows 

/■ / Vujx + z) \ ( Vu{x) \ ( Vu{x) \ Vu{x) 

"^'^' "^ " JuM.^.^^H.n v^'^ ■ \^n{x+z)\) ^^^^^ '^ [w^\) [w^^ij ■ w^m 

(5) 
where a denotes the surface measure on the "curve" {z : u{x + z) = u{x)} and where 
^G = /{^:v«(x)-z=o} G{z)a{dz) is a vector field of M^. 



Remark that the example we gave above is of this form. Indeed 



It is quite clear on this new formula that the singular integral operator is geometric 
(in the sense that it only depends on the curve and not and its parametrization u) and 
"fractional" . 

After this work was finished, we have been told that non-local minimal surfaces are 
being studied by Caffarelli, Roquejoffre and Savin [lO]. Loosely speaking, they study sets 
whose indicator functions minimize a fractional Sobolev norm || ■ ||//q, a G (0,1). They 
prove in particular that local minimizers are viscosity solutions of fi;[x,M] = 0. 

Comments and related works 

We gave two different formulations in the case of singular potentials. We think that 
Formulation (jlj) is the proper one in order to get a complete level set formulation of the 
geometric flow even if Formulation ([5]) is somehow more intuitive since it only involves 
the curve itself. In particular, the approach proposed by Slepcev [26] can be adapted (see 
dUD below). 

The level set equation we study has the following form 



dtu = fi{Du) 



Ci(x) + fi:[x,M] 



\Du\ in(0,+oo)xM^ (6) 



supplemented with the following initial condition 

u{0,x)=uoix) inM^ (7) 

where p denotes p/\p\ ii p ^ 0, n denotes the mobility vector field, ci(x) is a driving force. 
Equation ([H]) is a non-linear non-local Hamilton- Jacobi equation. A lot of papers are 
dedicated to the study of such equations. In our case, the main technical issues are the 
definition of viscosity solutions, the proof of their stability and the proof of a strong 
uniqueness result. We somehow use ideas from [26] and combine them with the ones from 
[7], even if the results of these two papers do not apply to our equation. 

From a physical point of view and as far as dislocation dynamics is concerned, the 
measure i/ldz) = co{z)dz should be uldz) = glz/lzDlzl'^^'^dz but in this case, the frac- 
tional mean curvature is not well defined (see Remark [T]). It is also physically relevant to 
say that close to the dislocation line, in the core of the dislocation, the potential should be 
regularized. On the other hand, it is important to assume that uldz) ~ g{z/\z\)\z\~^~^dz 
as \z\ -^ +CXD since this prescribes the long range interaction between dislocation lines. 
Another way to understand this difficulty is to say that in the core of the dislocation, the 
potential is very singular and the singularity should be compensated at a higher (second) 
order. On one hand, this can explain the loss of inclusion principle for such flows (if one 
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can define them for large times). On tlie otlier liand, one can tliink tliat in tliis case, 
tfie first term in sucli an expansion sliould be a mean curvature term. Tliis can make 
sense since curvature terms are commonly used to describe dislocation dynamics. It can 
be relevant to add one in ([6]). However, we choose not to do so in order to avoid techni- 
calities and keep clear some important points in the proof of the stability result and the 
comparison principle. 

In order to better understand properties of the fractional mean curvature flow, a 
deterministic zero-sum repeated game is constructed in [20] in the spirit of 



Organization of the article. In Section [21 we first give the precise assumptions we 
make on data. We next give the definition(s) of the fractional mean curvature k[x, ■]. In 
Section [31 we first give the definition of viscosity solutions for ([6]), we then state and prove 
stability results. We next obtain strong uniqueness results by establishing comparison 
principles. We also construct solutions of ^ by Perron's method. We finally give two 
convergence results which explain in which limit one recovers the classical mean curvature 
equation. In Section [H we verify that the zero levet set of the solution u we constructed 
in the previous section only depends on the zero level set of the initial condition. This 
provides a level set formulation of the geometric flow. In the last section, we give an 
alternative geometric definition of the flow in terms of generalized flows in the sense of 

Notation. S^^^ denotes the unit sphere of M^. The ball of radius 5 centered at x is 
denoted by Bs{x). If x = 0, we simply write Bs and if moreover 6 = 1, we write B. If 
p e M^ \ {0}, p denotes p/\p\. If A is a subset of M"' with d = N,N + 1 for instance, 
then A'^ denotes M'^ \ A. For two subsets A and B, AU B denotes AU B and means that 
A n i? = 0. The function 1a{z) equals 1 if z E A and if not. 

Acknowledgements. This paper is partially supported by the ANR grant "MICA". 
The author thanks R. Monneau and P. E. Souganidis for the fruitful discussions they had 
together. He also thanks G. Barles and the two referees for their attentive reading of this 
paper before its publication. 

2 Preliminaries 

In this section, we make precise the assumptions we need on data and we give several 
definitions of the fractional mean curvature. 

2.1 Assumptions 

Here are the assumptions we make on the singular measure throughout the paper. 

Assumptions. 

(Al) The mobility function /x : S^~^ -^ (0, +oo) is continuous. 



(A2) The driving force Ci : M^ ^ M is Lipschitz continuous. 

(A3) The singular measure i^ is a non-negative Radon measure satisfying 

' for all 5 > 0, u{R^ \ Bs) < +00 , 

for all r > 0, e e S^^\ iy{z e B : r\z ■ e\ < \z - {z ■ e)e\'^} < +cx) , 

6u{R^ \Bs)^0 as (5 ^ , 
for all e G S^~^, r v{z G B : r\z ■ e\ <\z — [z ■ e)ep} ^0 as r — »• 

(8) 
{B^ denotes the ball of radius 5 centered at the origin and B = Bi) and the last 
limit is uniform with respect to unit vectors e G S^~^. 

(A4) The initial datum uq '■ M^ — i> M is bounded and Lipschitz continuous. 

We point out that the set {z & B : r\z ■ e\ < \z — {z ■ e)ep} appearing in the second and 
the fourth lines of ([H]) is the region between an upward and downward (with respect to 
vector e) parabola. 

Even if the assumptions on the singular measure look technical at first glance, they 
are quite natural in the sense that they imply several important properties: 

• the measure is bounded away from the origin; 

• the singularity at the origin (if any) is a weak singularity in the sense that the 
fractional mean curvature of regular curves can be defined; if the reader thinks of 
the example given in ([3]), this means that we choose a < 1; 

• Parabolas {z : rz^- = {z'l"^} (which are the model regular curves for us) can be 
handled, even when they degenerate (r -^ 0). 

Example 1. The Standing Example for the singular measure is 

PSEidz) = g l^j —^ 

with g : S^~^ -^ (0, +00) continuous and a G (0, 1). The measure in ([3]) corresponds to 
the isotropic case {g = 1). 

2.2 Fractional mean curvature 

In this subsection, we make precise the definition of fractional mean curvature. Our 
definition extends the ones given in [121 [IZj where u^dz) = Co{z)dz to the case of singular 
measures. 

Let us define the fractional curvature of a smooth curve T = {x E M^ : u{x) = 0} = 
d{x G M^ : u{x) > 0} associated with u. If u is C^'^ and Du{x) ^ 0, then the following 
quantity is well defined (see Lemma [T] below) 

K*[X, F] = K*[X,U] = k\[x,u]- K-[X,U] ,^. 

K^[X,T] = K^[x,u] = K^[x,u] — K*_[x,u] 
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Q = {u > u{x)} 




{u < u{x)} 



T = {u = u{x)} 



where 



and 



Figure 1: Fractional mean curvature of a curve 



K^[x,u\ = u[^z : u{x + z) > u{x), Du{x) ■ z <Q) 
n~[x,u] = u[^z : u{x + z) < u{x), Du{x) ■ 2; > O) 



(10) 



K^[x, -u] = u^z : u{x + z) > u{x), Du{x) ■ 2 < O) 
K*_[x,u\ = u(^z : u{x + z) < u{x), Du{x) ■ z > O) . 

We will see later (see Lemma [H] below) that these functions are semi- continuous and this 
explains the choice of notation we made. In order to understand the way these quantities 
are related to the geometry of the curve {u = u{x)}, it is convenient to write for instance 

k'^[x, u] = v^z : < —Du{x) ■ z < u{x + z) — u{x) — Du{x) ■ z) . 

As shown on Figure [U, n^ [x, u] measures how concave the curve is "near" x and k~ [x, u] 
how convex it is. 

Lemma 1 (Fractional mean curvature is finite). If u is C^'^ at point x, i.e. there exists 
a constant C = C{x) > such that for all z G M^ 

\u{x + z) — u{x) — Du{x) ■ z\ < C\z\'^ 

and its gradient Du{x) 7^ 0, then k^[x, u] are finite. 

Ifu is C^'^ at X and Du 7^ everywhere on {y G M^ : u{y) = u{x)} and v is absolutely 
continuous with respect to the Lebesgue measure, then /t;^[x, u] are finite and 



K [x, u\ = K,^[x, u\ . 
Remark 1 . One can check that this lemma is false if a = 1 in the Standing Example [H 



Proof. We only prove the first part of tlie Lemma since the second part is clear. 
Since u is bounded on M^ \ Bg for all 6 > 0, it is enough to consider 



K*Y'^\X,U] = V 



V 



(^z E Bs : u{x + z) > u{x), Du{x) ■ z < O) 

(^z E Bs : < re ■ z < u{x + z) — u{x) + re ■ z) 



where r = |Z}M(a;)| ^ and e = r ^Du{x). If now zn denotes e ■ z and z' = z — z^e, and 
if we choose S such that r — C6 > 0, we can write 



K^) ' [X,U\ 



< iy{z e Bs:0<rzN <Czjf + C\z'\ 



< u{zeB5:0<C-\r-C6)zN<\z'\'^) 

and the result now follows from Condition ([8]). D 

The following lemma explains rigourously the link between (jH]) and ([2]) and the link 
with the formulation used in [TT] in the case where z/ is a bounded measure. 



Lemma 2 (Link with regular dislocation dynamics). Consider cq G L^(M^) such that 
Co{x) = Co{—x) . Then 

co{z)dz = - co + K^[x,u\ 

{z:u{t,x+z)>u{t,x)} ^ J 

sign* {u{x + z) — u{x))co{z)dz = -k*[x,u] 
sign^{u{x + z) — u{x))co{z)dz = -k*[x,m] 

with sign*(r) = 1 (resp. sign^(r) = 1) if r > (resp. r > 0) and —1 if not and with 
p{dz) = co{z)dz. 

Since the proof is elementary, we omit it. 

We conclude this section by stating two results which explain the link between two 
special cases of fractional mean curvature operator and the classical mean curvature op- 
erator. The first one appears in [1^] (see their Corollary 4.2). We state it in a special case 
in order to simplify the presentation. 

Proposition 1 (From dislocation dynamics to mean curvature flow - [IS]). Assume that 
V = v^ has the following form 



v{dz) = u%dz) = ^^+iM^^| Co (^) dz 



with Cq even, smooth, non-negative and such that Cq{z) = \z\ ^ ^ for \z\ > 1. 
Assume that u E ^^(M^) and Du{x) ^ 0. Then 

k[x,u\ = K^[x, m] -^ Cdiv I I (x) 

as e -^ for some constant C > 0. 
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Remark 2. In [To], general anisotropic mean curvature operators can be obtained by 
considering anisotropic measures u{dz). 

This result can be compared with the following one. 

Proposition 2 (From fractional mean curvature to mean curvature). Assume that v has 

the following form 

dz 
v{dz)=v-{dz) = {l-a)j^^ 

with a G (0, 1). 

Assume that u E C^iR^) and Du{x) ^ 0. Then 

k[x, m] = K°'[x,u\ — >• Cdiv 



as a -^ 1, a <1, where C is some positive constant. 

Remark 3. Anisotropic mean curvature can be obtained by considering 

."(d.) = (l-a)j(i|)i^. 

Sketch of the proof of Proposition^^ For all t], we first choose S such that 

\u{x + z) — u{x) — Du{x) ■ z D'^u{x)z ■ z\ < ri\z\'^ . (11) 

If e denotes —Du{x) and W{z) denotes u{x + z) — u{x) — Du{x) ■ z, we have 

fi:"[x,M] = u^'izeR^ ■.0<e-z<W{z)} 
-z/"{^GM^ : W{z) <e-z<0} 

dz . r dz 

A — a] / , ,,,, (1 — a 



{z(^Bs:Q<e-z<W{z)} 1^1 J {z(^Bs:W{z)<e-z<Q} Fl 

+0(1 -a) 

since |2;|~^~" is a bounded measure in B^. 

In view of (1111) . it is enough to prove the result for W{z) = Bz ■ z where i? is a 
symmetric N x N matrix. Hence we study the convergence of 

K'^ = (I -a) I -^-(1-a ^ "^^ 



{z£Bs:0<e-z<Bz-z} l^^l J {zeBs:Bz-z<e-z<0} Fl 

We next use the following system of coordinates: zi = e- z and z = {zi, z'). We now write 

Bz ■ z = bizf + Zi(b[ ■ z) + B' z ■ z 

for some hi G M, h'^ G R^^^ and a [N — 1) x [N — 1) symmetric matrix B' . We thus want 
to prove 

i^" -^ \e\^hvB' 

11 



as a ^ 1. We can assume without loss of generality that |e| = 1. For z G -85, we have 

e-z<Bz-z => zi<{l-C6)-^B'z' -z' 
^1 > (1 - C5)'^B'z' -z => e-z>Bz-z. 

Hence, it is enough to study the convergence of 

K'' = (l-a) I t4^ -il-a ^ "^^ 



{izi,z')eBs:0<zi<B'z'-z'} 1^1 "*"" J {(zi,z')eBs:B'z'-z'<zi<0} \^\ "*"" 

If a{d9) denotes the measure on the sphere S^~^, we can write 



j{(zi,z'):\z'\<5,0<zi<B'z'-z'} 1^1 "*"" J {(zi,z'y.\z'\<5,B'z'-z'<zi<0} \^\ "*"" 

Jg^SN-2.B'0.0>oJr=oJzi=Q 1^1 + '^ J^ ^' 

/■ /■'5 /-O ^Ar-2 

-{I -a) / / / 2 , „2NrAr^.V2 ^(^^)^^^^l • 



2+^2)(iV+a)/2' 



je&^-'^:B'e-e<o Jr=o J zi=r^B'e-e [^ 
We next make the change of variables zi = r'^r and we get 

Je&«-^:B'e-e>o Jr=o Jt=o [r t -\- i)'- ^ " 

We finally remark that 

V^ ^ (Q^ ^)^ £J ' (,2,2 + l)(^+.)/2 - i^^g ■ g as ^ ^ ■ 
In particular, for 6 small enough, 

(l-^)^^^-^<£f (,2,2 + l)(^+.)/2 <a + ^)^^^-^- 

It is now easy to conclude by remarking 

(1-a) f r-'^ = 5^-" 



'TV-l 



/" 9 9a {(19) = Ch 
where / denotes the (iV — 1) x (A^ — l)identity matrix and C is a positive constant. D 
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3 Viscosity solutions for ([6]) 

3.1 Definitions 

The viscosity solution theory introduced in [26] suggests that the good notion of so- 
lution for the fractional equation ([6]) is the following one. 

Definition 1 (Viscosity solutions for (E])). 1. An upper semi- continuous function u : 
[0, T] X R^ is a viscosity subsolution of ([6]) if for every smooth test-function such 
that u — (f) admits a global zero maximum at {t,x), we have 



dt<t){t,x) <fi{D(P{t,xj) 



Ci(x) + K*[x,(f)(t,-)] 



\Dcf>\{t,x) 



(12) 



if D({)[t, x) 7^ and dt4>{t, x) < if not. 



2. A lower semi- continuous function u is a viscosity supersolution of IQ if for every 
smooth test-function cj) such that u — cj) admits a global minimum at (t, x), we have 



dt(p{t,x) >fi{D(P{t,x)) 



Ci(x) + K^[x,(j){t,-)] 



W\{Xo,to) 



(13) 



if D(f)(t, x) ^ and dt4>(t, x) > if not. 

3. A locally bounded function u is a viscosity solution of (El) if u* (resp. u^) is a 
subsolution (resp. supersolution) . 

Remark 4. Given 5 > 0, the global extrema in Definition [1] can be assumed to be strict 
in a ball of radius 5 centered at (t, x). Such a result is classically expected and the reader 
can have a look, for instance, at the proof of the stability result in [7j. 

If one uses the notation introduced in [26], the equation reads 

dtu + F(x, Du, {z : u{x + z) > u{x)}) = (14) 

i^ and K C M^, 



with, for x,p E 
F{x,p,K) = 



-f^iP) 



ci{x) + v{K n {p ■ ^ < 0}) - ij{K'' n {p ■ ^ > 0} 







\p\ if p 7^ 
if not 



where K'^ is the complementary set of K. With this notation in hand, one can check that 
this non-linearity does not satisfy Assumption (F5) of [26]. The idea is to check that, 
somehow. Assumption (NLT) in [7] is satisfied and stability results thus hold true. 

Let us be more precise. We previously associated with n[-, ■] the following non-local 
operators (see the proof of Lemma [Tj) 



K 



+ \l,5\ 



,-+^2,5^ 



X,0] 



= v[z G Bs : (t){x + z) > (f){x), z ■ D(j){x) < O) 
[K-;p''[x,p,(j)] = u{z ^ Bs ■.(l){x + z) > (j){x), z-p<0) . 

In the same way, we can define 



(15) 
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the negative non-local curvature operators {k'Y'^, i = 1,2, 

upper semi-continuous envelopes of these four integral operators (k^)*'*^, i = 1,2 

and lower/upper semi-continuous total non-local curvature operators {k^, 



i,5 ( ^.*\i,5 



5 



1.2. 



By using the idea of Lemma [21 it is easy to see that 

K*f'^[x,p,u] = u(^z ^Bs ■.u{t,x + z) >u{t,x)'j -u{z ^Bs -.p- z > 0) , , . 



{k^:Y'^[x,p,u] = u[z ^ Bs : u{t,x + z) > u{t,x)) — u{z ^ Bs : p ■ z >0) 
We can now state an equivalent definition of viscosity solutions of ([6]). 

Definition 2 (Equivalent definition). 1. An upper semi- continuous function u : [0,T]x 
M^ is a viscosity subsolution of ([6]) if for every smooth test-function (p such that 
u — (f) admits a maximum at {t,x) on Bs{t,x), we have 



dt<t){t,x) <p{D(P{t,x)) 



C,{x) + {Ky''[x,<P{t,-)] + {Ky''[x,u{t,-)] 



\D<P\{t,x) (17) 



if D(j)(t, x) ^ and dt4>(t, x) < if not. 



2. A lower semi- continuous function u is a viscosity supersolution of ([6]) if for every 
smooth test-function (f) such that u — (p admits a global minimum at [t, x), we have 



dt(p{t,x) >^i{D(t){t,x)) 



Ci{x) + {Ky^'[x,(P{t,-)] + {Ky''[x,u{t,-)] 



\D<t>\{x^,t^) {li 



if D^{t, x) 7^ and dt4>(t, x) > if not. 



3. A continuous function u is a viscosity solution of i^ if it is both a sub and super- 
solution. 

Remark 5. Equivalent definitions of this type first appeared in [23] and since the proof is 
the same, we omit it. 

Remark 6. Remark H] applies to the equivalent definition too. 

Remark 7. Definition [2] seems to depend on 6. But since all these definitions are equivalent 
to Definition^ it does not depend on it. Hence, when proving that a function is a solution 
of (jH]), it is enough to do it for a fixed (or not) 6 > 0. 

3.2 Stability results 

Theorem 1 (Discontinuous stability). Assume (A1)-(A3). 

• Let {un)n>i be a family of sub solutions of ([6]) that is locally bounded, uniformly with 
respect to n. Then its relaxed upper limit u* is a subsolution of (EI). 
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• // moreover, u„,(0, x) = Uq{x), then for all x G M^ 

u*{0,x) <u*q{x) 
where u^ is the relaxed upper limit o/uq. 

• Let {ua)a<^A ^6 a family of suhsolutions of ^ that is locally hounded, uniformly 
with respect to a & A. Then u, the upper semicontinuous envelope of sup^, Ua is a 
subsolution of ([6]). 

Even if this result follows from ideas introduced in [7j together with classical ones, we 
give a detailed proof for the sake of completeness. 

Proof. We only prove the first part of the theorem since it is easy to adapt it to get a 
proof of the third part. The second one is very classical and can be adapted from [1] for 
instance. 

Consider a test function ip such that u* — ip attains a global maximum at (t,x). We 
can assume (see Remark Hj) that u* — ip attains a strict maximum at {t,x) on Bs{t,x). 
Consider a subsequence p = p{n) and (tp,Xp) such that 

u*{t,x)= lim Up/n)(tp,Xp) . 

n—*+co 

Classical arguments show that Up — ip attains a maximum on Bs(t, x) at {sp, yp) G Bsit, x) 
and that 

(■Sp, Vp) -* {t, x) and Up{sp, yp) -^ u*{t, x) . 

Since Up is a subsolution of ([6]), we have 



dt^pisp, yp) < 
fi{Dip{sp,yp)) 



ci{yp) + {K.*) ' [yp,(p{spr)] + {K.*) ' [yp,D^ip{sp,yp),u{spr)] 



\D<p\isp,yp) 



if Dip{tp, Xp) ^ and dt^pitp, Xp) < if not. If there exists a subsequence q oi p such that 
Dip{sq, yq) = 0, then it is easy to conclude. We thus now assume that Dip{sp, yp) ^ 
for p large enough. In view of the continuity of [i and ci, the following technical lemma 
permits to conclude. D 

Lemma 3. Assume that Dip{sp^yp) ^ forp large enough. 

• Assume moreover that D(p(t,x) ^ 0. Then 

is,y)^ {K*)^'^[y,ip{s,-)] and {s,y) h^ {K*)'^'^[y, D,^(p{s,y),Up{s, ■)] 

are well defined for i = 1,2 in a neighbourhood of (t, x) and 

]imsnpi{K*)'''[yp,^{sp,-)]\ < {K*)'''[x,^{t,-)] 

limsup I {K*)'^'^[yp, , D^(p{sp, yp), u{sp, ■)] i < {k*)'^'^[x, D^ip{t, x), ip{t, ■)] 
as soon as Up{sp, yp) -^ u{t, x) as p —^ +oo. 
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Assume now that Dip{t, x) = 0. Then, for i = 1,2, 



f^*) ' [yp,^{sp,-)] + {K-*) ' [yp,Dip{sp,yp),u{spr)] 



\Dip\{sp,yp) -^0 asp^+oo. 



As we shall see, this lemma is a consequence of the following one. 

Lemma 4 ([26]). Consider fp and Qp two sequences of measurable functions on a set U 
and f > limsup* fp, g > limsup* gp, and ap, bp two sequences of real numbers converging 
to 0. Then 

^{{fp > ap, Qp > ^p} \ {/ > 0, 5- > 0}) ^ as n^ +oo . 

We mention that in [26], the measure is not singular and there is only one sequence of 
measurable functions but the reader can check that the slightly more general version we 
gave here can be proven with exactly the same arguments. An immediate consequence of 
the lemma is the following inequality 

lim sup z/({/p > ap,gp > bp}) < z/({/ > 0,5- > 0}) . 

p 

Proof of Lemma\^ Let us first assume that Dif(t,x) ^ 0. In this case, for {s,y) close to 
(t,x), Dif{s,y) 7^ and all the integral operators we consider here are well defined (see 
Lemma d]). Recall next that, for i = 1,2, (/t*)*''' = (k^^)*''^ — («:")*'''. Hence, it is enough 
to prove that 

limsup |«)i''^[l/p,y.(sp,-)]| < {^lY''[x,^{t,-)], 

\immih^:Y^'[yp,^{Sp,-)]\ > {^:Y''[x,^{t,-)] , 

limsup <^ {^If'^lyp, D^ip{sp,yp),Up{sp, ■)] \ < {K\f'^[x,D^(f{t,x),u*{t,-)], 

liminf <^ {K~Y'^[yp,D^(f{sp,yp),Up{sp,-)] > < {n-f'^lx, D^(f{t,x),u*{t, ■)] . 

In order to prove the first inequality above for instance, choose fp{z) = ip{sp, yp + z) — 
(p(t,x), Qp = (p{sp,yp) - ip{t,x), gp{z) = -D(p{sp,yp) ■ z, bp = in Lemma H 

We now turn to the case Dipit, x) = 0. We look for 6 = 6p that goes to as p ^ +oo 
such that \Dip{sp,yp)\ < CSp and 

{Kiy'^''[yp,ip{sp,-)]\D^{sp,yp)\ ^0 and {K;y'^-[yp,ip{spr)]\D^{sp,yp)\ ^0 

as p ^ +00. This is enough to conclude since Condition ([H]) implies that 

{Klf'^^lyp, Dip{sp,yp),u{sp, ■)]\Dip{sp,yp)\ ^ 
(i^*)'^'^'' [Vp, Dip{sp, yp),u{sp, ■)]\Dip{sp, yp)\ ^ . 
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We only prove that the first hmit equals zero since the argument is similar for the second 
one. If Tp denotes \Dip[sp,yp)\ and Cp denotes —r~^Dip{sp,yp), and z^ = ep ■ z and 
z' = z — z^Cp, then 

{KlY'^lyp, ip{sp, ■)]\Dip{sp,yp)\ 

= rpu{z e Bsj, -.0 < TpCp ■ z < ip{sp, yp + z) - ip{sp, yp) + VpCp ■ z) 
< rpv{z eBs^:0< VpZN < C\z'\^ + Cz^) 

< rpv{z e Bsp :0 < VpZ^ < C\z\^ + CbpZj^) 

where C is a bound for second derivatives of c/? around (t, x). Now if we choose 6p = 

rp/{2C), we get 

{K\Y^'[yp,ipiSpr)]\Dipisp,yp)\ < Vpuiz E Bs^ : < irp/2C)z^ < \z'\^) 

< rpu{z G 5 : < {rp/2C)zN < |;zT) 

and the last limit in ([8]) permits now to conclude. D 

3.3 Existence and uniqueness results 

Let us first state a strong uniqueness result. 
Theorem 2 (Comparison principle). Assume (A1)-(A4)- Assume moreover 
(A3') For all e e S^"^ and r e (0, 1) 

r u{z e Bs : r\z ■ e\ < \z - {z ■ e)e\'^} ^ as 5 ^ (19) 

uniformly in e and r G (0, 1) and 

v{dz) = J{z)dz with J G iy^'^(M^ \ Bs) for all 6 > . (20) 

Consider a bounded and Lipschitz continuous function uq- Let u (resp. v) he a hounded 
suhsolution (resp. hounded supers olution) of ([6]). Ifu{0,x) < uo{x) < v{0,x), then u <v 
on (0,+cx)) X M^. 

The proof is quite classical. The main difficulty is to deal with the singularity of the 
measure. 

Proof of Theorem\^ We classically consider M = sup^ ,j,{M(t,x) — v(t,x)} and argue by 
contradiction by assuming M > 0. We next consider the following approximation of M 

Ms^a = sup {u[t,xj — v{s,yj e — rjt — a\x\ }. 



t,s>0,x 



27 2e 
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Since u and v are bounded, this supremum is attained at a point {i,s,x,y). We first 
observe that M^^a > M/2 > for 77 and a small enough. Since u and v are bounded, this 
implies in particular 

r/t + e^*^^ + a|x|2<Co (21) 

where Cq = IklU + Halloo- 

Classical results about penalization imply that (t, s, x, y) -^ (t, t, x, |/) as 7 ^ and 
(t, i,x,y) realizes the following supremum 

I 12 
M^^a = sup {u{t,x) — v{t,y) — e^^ rjt — a\x\'^} . 



It is also classical [M] to get, for e and 77 fixed, 

a|x|2 ^0 asa^O. (22) 

We claim next that this supremum cannot be achieved at t = if e, a,?7 are small 
enough. To see this, remark first that M^^a > M/2 > for rj and a small enough and, if 
t = 0, use the fact that uq is Lipschitz continuous and get 

J\/J 1^ — 7/1 T 1 

< — < sup {uoix) - uoiy) } < supjCor - — } = -CqE 

■^ x,y&RN ^^ '-■>0 '^^ ^ 

and this is obviously false if e is small enough. We conclude that, if the four parameters 
are small enough, t > and s > 0. 

Hence, we can write two viscosity inequalities. In order to clarify computations, we 
introduce the function M(p) defined as follows 

^^^ \ if |) = . 

It is easy to see that M is uniformly continuous and it trivially satisfies 

\M{v)\< ||/i||oob|. 

In the following, ujm denotes the modulus of continuity of M. 
We now write viscosity inequalities: for all 5 > 0, 

r/ + ^^ + i^e^*^^i < (ci(x) + {K*f'\i, Ml ■)] + («:*)''^[x,p + 2ax, u{i, ■)])M{p + 2ax) 

7 Ze 



7 
where p = e^*^^ and 



^ ' > {cM + {^.Y''[y,Us,-)] + {^,)'''[y,p,vCs,-)])Mip) 



(j)u{t,x) = v{s,y) + ^^-^ + e^'^-^^-—^ + r]t + a\x\^, 

27 2e 

_ [ e t) ~\v x\ ~ 

(t)v{s,y) = u{t,x) e^* — r]t-a\x\^. 

27 Is 
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Substracting these inequalities yield 

I ~ ~|2 I ~ ~|2 

n + Ke^^^Zl^ < ll^ll^pdlUe^*^^^ + \\c,\\^ujMi2^/C^) + T^i (23) 

le £ 

(we used (ED) where 

(«:,)i'^[y,0,(S, ■)] + {^^.f^\y.V.v{s, ■)]]m{P) . 



Our task is now to find 5 = 6{a,e) so that the right hand side of this inequality is small 
when the four parameters are small. We distinguish two cases. 

Assume first that there exists a sequence a„ ^ and £„ — > such that 

p = pn^O. 

In this case, we simply choose S = 1, K = 2||/i||oo||-Dci||oo and we pass to the limit as 
n —* +00 in ( |23|) and we get the desired contradiction: f] < 0. 

Assume now that for a and e small enough, we have a constant C^ independent of a 
such that 

|p|>C, >0. (24) 

In this case, the following technical lemma holds true. 

Lemma 5. By using { \19\i . we have 

Tni < -os{l) + La/(2v/CW) + o^{l)[e] + C,e^*^^^ 
e e 

where Cq appears in fl^T]) and Cs only depends on J, ||/i||oo one? 6 (we emphasize that the 
third term goes to as a -^ for fixed e). 

The proof of this lemma is postponed. We thus get (recall that p = e^*(x — y)/£) 

rj+Ke^it^ < Ce^*-^^^+C(l + l)c.M(2v/CW) + -05(l)+o.(l)[e]+C,e^*^^^ 
2e e e e 

where C only depends on Ci, v and ||m||oo + ll'^Hoo and Cs is given by the lemma. By 
choosing K = 2{C + C5), we get 

V<C{1 + \)ujm{2^C^) + -0,(1) + o^{l)[s] . 
e 

By letting successively a and 5 go to 0, we thus get a contradiction. This achieves the 
proof of the comparison principle. D 
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Proof of Lemma O We first write 

Tni < \\fi\\oo\{fi*Y'^\[x, (puii, ■)]]? + '^ax\ + \\fi\\oo\{fi*y'^\[y,(f)v{s,- 
+ \iK*f'^\[S:,P + 2ax,uii,-)]ujMi\2ax\) 

+ ({K*y''[x,p + 2ax, u{i, ■)] - {n,y^'[y,P, v{~s, ■)]] M{p) . 



We thus estimate the right hand side of the previous inequality. We start with the first 
two integral terms. 

\{K*y^'\[x,Ui,-)] < KY''[i,Ui,-)] + KY^'[x,Ui,-)] 

< u{ze Bs:0< -{p + 2ax) ■ z < {a + 6^7(2^))^^) 
+u{ze Bs:0> -{p + 2ax) ■ z > {a + e^^^/{2e))\z\'^) 

< iy{z G Bs : \ep + 2eax\\e- z\ < C{r])\z\'^) 

where we use (^T^ to ensure, for a,e small enough, 

Ki 1 

— + «£<-e^''"/^ + l:=C(r/). 
If now Ta^e deuotes \ep + 2eax\ and we choose 6 < r ^^^ / {2C {rj)) , we use flT^ to write 

|(K*)^'^|[x,0„(t,-)]|p + 2ax| = -rc,^eT^{zeBs:r^^e\e-z\<C{vi)\z\^) 

< -ra,ei^{z e Bs : -ra,e\e ■ z\ < C{r])\z - (e ■ z)e\^) 

2C(r/) r . ^ , , , . X r 

< < sup rz/(z G Bs : r\e ■ z\ < \z — [e ■ z)e\ 

= ^osil). 

Since ax — > (see (^^), we choose for instance 

e\p\ 



6< 



4C(r?) 
Arguing similarly, we get for 6 < -^K, 

\{K^Y^'\[y,MSr)]\P\<losil). 
As far as the third integral term is concerned, we simply write 

\iK*y^'\[x,uii,-)]uMi\2ax\) < iy{BI)ujMi\2ax\) < liJMi2^/C^) 



20 



(we used f l^T]) ). We now turn to the last two integral terms. In view of ( 1161) . we can write 

f,z = {K*f^'[x,p + 2ax,u{i,-)]-{K,f^'[y,p,v{s,-)] 

= ^{z ^ Bs : u{i, X + z) > u{i, x)) — i/{z ^ Bs : v{s, y + z) > v{s, y)) 
-u{z ^ Bs : {p + 2ax) ■ z > 0) + i/{z ^ Bs : p ■ z > 0) . 

Now, we use (!2U|) to get 

Tnl = / J{Z - 5)l{„(t,)>n(t,x)}(^)c?^ - / y(^ - y)'i-{v{s,-)>u{s,y)}i^)dz + Oa{l)[e] ■ 

Remark next that the definition of (t, S, x, y) implies the following inequality: for all 

I ~ ~I2 

u{t, z) — u{t, x) < v{s, z) — v{s, y) + a{\z\^ — \x\^) — e^* — . 

This implies that for \z\ < Ra,ei we have 

l{«(t,-)>«(t-S)}(^) - '^{v{s,-)>u{s,y)}{z) 

where 

a \ Is 



> 






e 



where C^ appears in ( !24|) . We used here ( l22l) . Hence, we have 



Tnl < / J(2 — x)(i2 + / \J{z ~ x) — J{z — y)\dz 

< / J{z)dz + Cs\x — y\ = Oa{l)[£] + Cs\x — y\ 

J\z\> ^^^ 



■2^80(77)" 

where we used once again ( l22i) . It is now easy to conclude. D 

We now turn to the existence result. 

Theorem 3 (Existence). Assume (A1)-(A4) and (A3'). There then exists a unique 
bounded uniformly continuous viscosity solution u of ^, ([7]). 
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Proof. We first construct a solution for regular initial data. Precisely, we first assume 
that Mo G C^{M.^) (the function and its first and second derivatives are bounded). 

Because we can apply Perron's method, it is enough to construct a sub- and a su- 
persolution u"^ to ([6]) such that (m"'")*(0,x) = (m~)*(0,x) = uo{x). We assert that 
M^(t,x) = uo{x) ± Ct are respectively a super- and a subsolution of ([6]) for C large 
enough. To see this, we first prove that there exists Cq = Co(||-D^Mo||oo) such that for all 
X G M^ such that Duq{x) 7^ 0, we have 

{\k,\[x,uo] + \k*\[x,uo])\Duo{x)\ < Co. (25) 

In order to prove this estimate, we simply write for x such that Duq{x) 7^ 

((fi:;)i'^[x,Mo] + {<_f'^[x,UQ])\Duo{x)\ < 2v{z G Bs : r\e ■ z\ < ]^\\D^u4^\z\'')r 

< 2u{z G Bs : r\e ■ z\ < C\z - (e ■ z)e\'^)r < C^ 

where r = \Duo{x)\, C = max(||DMo||oo, 1) and e = Duo{x)/r and 6 = r/{2C) and Cy is 
given by (IHl). On the other hand 

(«)2'^[x,uo] + {K*_f^'[xM)\Duo{x)\ <^r = ICyC. 



We thus get Estimate ( l25l) . 

If now Mo is not regular, we approximate it with m'q G C|(]R^) and can prove that the 
corresponding sequence of solutions u^ converges locally uniformly towards a solution u. 
Since this is very classical, we omit details (see for instance PP). D 

We now explain in which limit one recovers the mean curvature fiow. To do so, we 
state two convergence results. Their proofs rely on Propositions [1] and El The first one 
(Theorem Hj) appears in [15] and the second one can be proved by using Proposition [21 

Theorem 4 ([IS])- Assume that fi = 1, ci = 0, Uq is Lipschitz continuous and bounded 
and 

iy{dz) = u'idz) = -^^^3^-^co (^) dz 

with Co even, smooth, non-negative and such that Cq{z) = \z\~^~^ for \z\ > 1. Then the 
viscosity solution m^ 0/ (Ej), ([7]) converges locally uniformly as e — ^ towards the viscosity 
solution u of 

dtu = C\Du\ div 



\Du 
(C is a positive constant) supplemented with the initial condition ((Tj) . 

Theorem 5. Assume that /x = 1, ci = 0, uq is Lipschitz continuous and hounded and 

dz 
u{dz) = u-{dz) = {l-a)—^ 
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with a G (0, 1). Then the viscosity solution n" of (^, ([7]) converges locally uniformly as 
a ^ 1 towards the viscosity solution u of 

dfU = C\Du\ div ' 



\Du\ 
(C is a positive constant) supplemented with the initial condition (I7j) . 

4 The level set approach 

In the previous section, we constructed a unique solution of ([6]) in the case of singular 
measures satisfying (A3) and (A3') and for bounded Lipschitz continuous initial data (see 
(A4)). In the present section, we explain how to define a geometric flow by using these 
solutions of ^. Precisely, we first prove (Theorem [H]) that if u and v are solutions of 
([6]) associated with two different initial data mq and vq that have the same zero level sets, 
then so have u and v. Hence, the geometric flows is obtained by considering the zero level 
sets of the solution u of ([6]) for any (Lipschitz continuous) initial datum. We also describe 
(Theorem [7]) the maximal and minimal discontinuous solutions of (E]) associated with an 
important class of discontinuous initial data. 

Theorem 6 (Consistency of the definition). Assume (A1)-(A3) and (A3'). Let uq and 
Vq be two bounded Lipschitz continuous functions and consider the viscosity solutions u, 
V associated with these initial conditions. If 



then, for all time t > 0, 



{x e M^ : Uq{x) > 0} = 
{a; e M^ : Uo{x) < 0} = 


= {xeR'' 


: Vo{x) > 0} 
: Vo{x) < 0} 


>0, 






{xeR^ ■.u{t,x)>0} = 
{xeR^ ■.u{t,x) <0} = 


= {xGM"^ 


■.v{t,x) > 0} 
: v{t,x) < 0} 



In view of the techniques used to prove the consistency of the definition of local 
geometric fronts (see for instance [8]), it is clear that this result is a straightforward 
consequence of the following proposition. 

Proposition 3 (Equation ([6]) is geometric). Consider u : [0, +cxo) x M^ a bounded sub- 
solution of (El) and 9 : R ^ R a upper semi- continuous non- decreasing function. Then 
9{u) is also a subsolution of (E])- 

Such a proposition is classical by now. It is proved by regularizing 9 (in a proper way) 
with a strictly increasing function ^", by remarking that K*[a;,^"(M)] = fi;*[x, m] in this 
case, and by using discontinuous stability. Details are left to the reader. 
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Thanks to Theorem El we can define a geometric fiow in the following way. Given 
(Fq, Dq, Dq) such that Fq is closed, D^ are open and R^ = Fq U D^ U Dq , we can write 

D^ = {x eR^ : Uo{x) > 0}, Dq ={xeR^ : Uo{x) < 0}, Fq = {x G M^ : Uo{x) = 0} 

for some bounded Lipschitz continuous function uq (for instance the signed distance func- 
tion). If u is the solution of ([6]) submitted to the initial condition u{0,x) = uo{x) for 
X G M^, then Theorem [6] precisely says that the sets 

D+ = {x G M^ : u{t, x) > 0}, D- = {x G M^ : u{t, x) < 0}, F* = {x G M^ : u{t, x) = 0} 

does not depend on the choice of uq. 

The next theorem claims that there exists a maximal subsolution minimal supersolu- 
tion of (jH]) associated with the apropriate discontinuous initial data. 

Theorem 7 (Maximal subsolution and minimal supersolution). Assume (Al)-(A3) and 
(A3'). Then the function l^i+^p — Ijj- (resp. 1^+ — '^j^-yj^ ) is the maximal subsolution 
(resp. minimal supersolution) of ([6]) submitted to the initial datum '^o+^jy^ ~ Id- (resp. 

This result is a consequence of Proposition [3] together with discontinuous stability and 
the comparison principle. See [8], p. 445] for details. 

We conclude this section by showing that a bounded front propagates with finite speed. 

Proposition 4 (Evolution of bounded sets). Assume (Al)-(A3) and (A3'). Let VLq be 
a bounded open set of M.^ : there exists R > such that Qq C Br. Then the level set 
evolution (Tt, Df, D^) of {OVLq, VLq, (^o)'^) satisfies Df U F^ C Br+ci with 

C=||ci||oo- inf z/(2 gM^ : < e-2 < |2p) 

gg§iV-l 

as long as R + Ct > 0. 

Remark 8. Another consequence of this proposition is that, if there are no driving force 
(ci = 0), then the set shrinks till it disappears. 

Proof. The proof consists in constructing a supersolution of (El), ((71). It is easy to check 
that C is chosen such that 



u{t, x) = Ct + Ve^ + i?2 - yi^ 



X 



2 



is a supersolution of (El). Since Br = {x G M^ : m(0,x) > 0}, we conclude that Df' U 
Fi C {x G M^ : u{t,x) > 0} = BR.^^t) with i?=(t) = ^{Ct + Ve^TW)'^ - e^. Hence, 
-Di U F( C r\e>oBR^{t) = Brj^ci- □ 
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5 Generalized flows 



In this section, we follow [9] and give an equivalent definition of the flow by, freely 
speaking, replacing smooth test functions with smooth test fronts. 

In order to give this equivalent definition, we use the geometrical non-linearities we 
partially introduced in Section [21 For all x,p E M.^ and all closed set JF c M^ and open 

set C C M^ 



F*{x,p,0) 



-KP) 



-f^iP) 



ci(x) + u^j" n{p-z<o})- ij{T'' n {j9 ■ ^ > 0} 
ci(x) + v{Or\{p-z < 0}) - v{0^ n {p ■ 2 > 0} 



\p\ if p 7^0, 
if not , 

IpI if p 7^ 
if not . 



We can now give the definition of a generalized flow. 

Definition 3 (Generalized flows). The family {Ot)te{o.T) of open subsets of M'^ (resp. 
{J^t)t&{o,T) of closed subsets ofMJ^) is a generalized super-flow (resp. sub-flowj of ([6]) if for 
all (to, xq) G (0, +cxd) X M^, r > 0, /i > 0, and for all smooth function : (0; +oo) x M^ -^ 
M such that 



1. dt(f) + F*{x, D(t), {z : 0(t, x + z)> (f){t, x)}) < -6^ in [to, to + h] x B{xq, r) 
(resp. dt(l) + F^{x,D(f), {z : 0(t,x + z) > 0(t, x)}) > -6^ in [to, to + h] x B{xo,r)) 

2. D(t)^0 m {(s, y) G [to, to + h]x B{xo, r) : 0(s, y) = 0}, 

3. {yGM^:0(to,y)>O}cO,\, 

(resp. {y G M^ : 0(to, y)<0}cR^\ J'J, 

4. {y ^ B{xo,r) : (f){s,y) > 0} C Ol for all s e [to,to + h], 

(resp. {y ^ B{xo, r) : 0(s, y) < 0} C M^ \ ^, for all s G [to, to + h]), 

then {y G B{xo,r) : 0(to + h,y) > 0} C Ol^^^ (resp. {y G B{xo,r) : 0(to + h,y) <0} C 

^ \ ^to+h)- 

Loosely speaking about generalized super-flows. Condition [1] says that in a prescribed 
neighbourhood V around (to,Xo), the normal velocity of the test front {0 > 0} is strictly 
smaller than the one of the front 0\ Condition [2] asserts that the front {</> = 0} is smooth 
in V; Conditions [3] and H] assert that the test front is inside the front O outside V. The 
conclusion is that the test front is inside the neighbourhood O at time t + h. 

Remark 9. As far as local geometric fronts are concerned. Conditions [3] and H] impose 
that the test front is inside O on the parabolic boundary of the neighbourhood. Here, 
because the front is not local, the test front has to be inside O everywhere outside the 
neighbourhood. 
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The next theorem asserts that Definition [3] of the fiow coincides with the level set 
formulation of Section |H 

Theorem 8 (Generalized fiows and level set approach). Assume (A1)-(A3) and (A3'). 
Let {Ot)t£(o,T) be a family of open subsets o/M^ (resp. {J^t)t(^{o,T) of closed subsets ofR^) 
such that the set Uig(o,T){0 ^ ^t is open in [0,T] x M^ (resp. Utfz{o,T){t} x JF^ is closed 
m [0, T] X M^;. 

Then {Ot)te{o,T) (resp. {J-'t)te{o.T)) is a generalized super-flow (resp. sub-fiow) of (El) if 
and only ifxi^^x) = lc)^(x) - lKiv\c,^(x) (resp. xi^^x) = 1jf^(x) - lKiv\^^(x)j is a viscosity 
supersolution (resp. subsolution) of ^, ([7]). 

Since the proof of [9] can be readily adapted, we omit it. We give a straightforward 
corollary of Theorems [7] and [H] that is used in 



Corollary 1 (Abstract method). Assume (Al)-(A3) and (A3'). Assume that {Ot)t (irid 
{J^t)t ore respectively a generalized super-flow and generalized sub-flow and suppose there 
exists two open sets Dq,Dq such that M^ = dOo U Dq U Dq and such that Dq C Oq 
and Dq C J-'q. Then if (F^, Df, D~[) denotes the level set evolution of ((90o, Dq , Dq), we 
have for all time t > 

Dt dOtd Dt U Vu D- C J'^ C D; C Vt . 

Remark 10. One can check that under the assumptions of the previous corollary, we have 
in fact Dq = Oq and -Dq" = JFq. 
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